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ABSTRACT 


The  United  States  Navy  is  at  a  crossroads  in  the  design  of 
ship's  engineering  plants.  Advances  in  solid-state  power 
electronics  combined  with  a  shift  to  gas  turbine  powered  propulsion 
and  electric  plants  has  placed  renewed  emphasis  on  developing 
advanced  power  systems. 

These  advanced  power  systems  may  combine  the  prime  movers 
associated  with  propulsion  and  electric  power  generation  into  an 
integrated  system.  The  development  of  advanced  electric 
distribution  systems  and  propulsion  derived  ships  service  (PDSS) 
power  systems  are  interim  steps  toward  the  goal  of  an  integrated 
system. 

Advances  in  the  design  of  ships  power  systems,  whether 
revolutionary  or  evolutionary,  will  require  extensive  testing  and 
simulation.  This  thesis  will  develop  a  basis  with  which  to  judge 
various  simulation  tools,  it  will  then  evaluate  a  simulation 
program  developed  for  the  Navy  by  the  Massachusetts  Institute  of 
Technology. 
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I .  INTRODUCTION 


A.  BACKGROUND 

The  United  States  Navy  is  at  a  crossroads  in  the  design  of 
ship's  engineering  plants.  Advances  in  solid-state  power 
electronics  combined  with  a  shift  to  gas  turbine  powered 
propulsion  and  electric  plants  has  placed  renewed  emphasis  on 
developing  advanced  power  system  3. 

These  advanced  power  systems  may  combine  the  prime  movers 
associated  with  propulsion  and  electric  power  generation  into 
an  integrated  system.  The  development  of  advanced  electric 
distribution  systems  and  propulsion  derived  ships  service 
(PDSS)  power  systems  are  interim  steps  towards  this  goal.  The 
end  result  may  be  an  Integrated  Electric  Drive  (IED)  which  is 
characterized  by  a  central  electric  generating  plant  providing 
electricity  both  to  synchronous  propulsion  motors  and  for 
ships  service  use.  The  driving  factors  in  this  program  are 
reduced  lifetime  costs,  weight  savings,  and  ship  quieting. 

B .  OBJECTIVES 

Advances  in  the  design  of  ship's  power  systems,  whether 
revolutionary  or  evolutionary,  will  require  extensive  testing 
and  simulation.  The  required  simulations  must  be  capable  of 
modeling  both  high  speed  electrical  transients  and  relatively 
low  speed  mechanical  transients.  Further,  they  must  be  able 
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to  effectively  model  fast  switching  power  converters  that  are 
an  essential  part  of  most  advanced  power  systems.  This  thesis 
will  develop  a  basis  with  which  to  judge  various  simulation 
tools.  Specifically  it  will: 

•  define  a  baseline  system  to  be  simulated  for  use  in 
judging  the  efficacy  of  various  programs 

•  list  a  set  of  evaluation  metrics 

•  detail  the  minimum  operating  conditions  to  be  modelled 

•  test  a  simulation  program  developed  for  David  Taylor 
Research  Center  (DTRC) 

C.  SUMMARY  OF  THESIS 

1.  A  View  of  Present  and  Future  Naval  Power  Systems. 

This  section  presents  the  basics  of  current  naval 
power  systems  and  will  detail  the  evolutionary  steps  currently 
envisioned  by  Naval  Sea  Systems  Command  Code  05Z  (NAVSEA  05Z)  , 
the  Advanced  Ship's  Machinery  Systems  Project  Office,  and 
David  Taylor  Research  Center  towards  the  goal  of  an  IED 
system. 

2.  Required  Software  Capabilities. 

The  required  capabilities  of  a  software  package  are 
compiled  as  specified  in  a  variety  of  sources.  Desirable 
capabilities  are  also  discussed. 
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3 .  Detailed  Models  Necessary  to  Simulate  a  Ships  Power 

System. 

This  chapter  provides  detail  on  specific  models 
developed  for  the  simulation  of  ships  power  systems. 
Mathematical  models,  block  diagrams,  or  schematics  are 
provided . 

4 .  Detailed  Description  of  WAVESIM 

The  computer  program  WAVESIM  developed  at  the 
Massachusetts  Institute  of  Technology  (MIT)  by  LT  Norbert 
Doerry  is  discussed  [1]  .  The  components  of  the  simulation 
program  are  described  and  the  solution  methods  used  by  Doerry 
to  address  the  problems  of  modelling  marine  power  systems  are 
summarized. 

5.  Analysis  and  Testing  of  WAVESIM 

This  chapter  presents  the  results  of  a  number  of 
simulations  run  using  WAVESIM.  It  then  presents  the  results  of 
an  evaluation  of  WAVESIM  using  the  metrics  developed  in 
Chapter  III.  This  chapter  is  not  a  technical  validation  of 
WAVESIM. 

6 .  Conclusions 

This  chapter  makes  a  statement  concerning  the  need  for 
a  program  suitable  for  simulating  advanced  marine  power 
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systems.  It  then  identifies  future  work  required  if  WAVESIM 
is  to  fill  that  need.  Lastly  it  identifies  other  simulation 
packages  that  might  be  useful  for  this  type  of  study. 
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II.  A  VIEW  OF  TEE  PRESENT  AND  FUTURE  OF  NAVAL  POWER  SYSTEMS 


The  current  design  of  naval  power  plants  is  similar  to 
that  used  for  many  decades.  Whether  a  ship  is  powered  by  steam 
turbines,  gas  turbines,  or  diesel  engines  the  prime  movers 
associated  with  propulsion  and  power  generation  are  separate. 
Relatively  unsophisticated  local  control  is  used  to  operate 
the  electric  power  plant.  Recent  advances  in  power 
electronics  and  control  of  power  systems  may  allow  significant 
advances  to  be  made  in  this  area  culminating  in  an  IED  System 
with  an  advanced  distribution  system. 

A.  PRESENT  DESIGN  OF  NAVAL  POWER  PLANTS 

Design  of  naval  power  plants  emphasizes  survivability  and 
weight  reduction.  Ships  power  distribution  systems  are  three 
phase  and  ungrounded  operating  at  a  nominal  450  volts. 
Electricity  is  typically  generated  by  two  or  three  synchronous 
generators  that  may  operate  in  parallel  or  isosynchronously . 
On  conventional  (non  nuclear)  ships,  the  generators  may  be 
driven  by  steam  turbines,  gas  turbine  engines,  or  diesel 
engines.  Table  I  lists  the  electrical  plant  characteristics 
of  a  selection  of  conventionally  powered  US  combatants. 
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TABLE  I 

ELECTRIC  PLANT  CHARACTERISTICS  OF 
SELECTED  US  SHIPS  [l:p.32] 

Ship  Class 

Class  Name 

Generation 

Prime 

Movers 

CG  16 

Lehey 

4  X  1500  KW 

1  X  750  KW 

1  X  300  KW 

Steam 

Turbine 

Gas 

Turbine 

Diesel 

FF  1052 

Knox 

3  X  750  KW 

1  X  750  KW 

Steam 

Turbine 

Diesel 

FFG  7 

Oliver  Hazard 
Perry 

4  X  1000  KW 

Diesel 

DD  963 

Spruance 

4  X  2000  KW 

Gas 

Turbine 

CG  47 

Ticonderoga 

3  X  2500  KW 

Gas 

Turbine 

DDG  51 

Arleigh  Burke 

3  X  3000  KW 

Gas 

Turbine 

Loads  are  categorized  as  either  vital  or  non-vital.  Vital 
equipment  has  at  least  one,  possibly  two,  redundant  sources  of 
power.  Equipment  essential  to  ships  safety  is  considered 
vital  including: 

•  Steering  motors 

•  Auxiliary  equipment  supporting  propulsion  and  power 
generation 

•  Damage  control  equipment  such  as  fire  pumps  and  interior 
communications 

•  Lighting 

•  Communications  and  navigation  equipment 

•  Weapons  systems 
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•  Machinery  space  ventilation 


Switching  of  vital  equipment  to  alternate  power  sources  is 
conducted  either  through  automatic  bus  transfer  (ABT)  switches 
or  manual  bus  transfer  (MBT)  switches.  Further,  the  type  of 
controller  found  on  vital  equipment  determines  whether  a  given 
piece  of  equipment  will  restart  automatically  when  power  is 
restored  as  is  the  case  for  a  low  voltage  release  (LVR) 
controller,  or  must  be  manually  restarted  as  is  the  case  for 
a  low  voltage  protection  (LVP)  controller.  Figure  1  shows  the 
generator  and  switchgear  configuration  for  a  typical  steam 
propulsion  ship  [2] .  Figure  2  shows  the  generator  and 
switchgear  configuration  for  a  modern  gas  turbine  powered  ship 
[3]  . 

B.  ANTICIPATED  DEVELOPMENTS  IN  SHIPS  POWER  SYSTEMS 
1.  Advanced  distribution  systems 

As  a  first  step  towards  an  I  ED  system  a  computer 
controlled  advanced  distribution  system  is  envisioned.  This 
system  would  have  the  ultimate  goal  of  automating  "monitoring 
and  control  decision  making  on  the  naval  power  systems  to  the 
greatest  extent  possible  [4]".  This  would  be  accomplished 
using  state  of  the  art  microcomputer  control  combined  with 
power  electronics.  Ideally  the  system  developed  will  have  the 
following  characteristics  [5]  : 

•  solid-state  bus  transfer 
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Figure  1  Generator  and  switchgear  arrangement  for  a  steam 
plant  [2] 
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Figure  2  Generator  and  switchgear  arrangement  for  a  gas 
turbine  plant  [3] 
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•  solid-state  circuit  breakers 


•  energy  storage 

•  reconfiguration  techniques 

•  embedded  intelligence 

A  number  of  the  points  listed  above  need  to  be  expanded  upon, 
a .  Energy  Storage 

The  number  one  design  issue  for  a  combatant  power 
system  is  continuity.  A  ship  must  maintain  its  ability  to 
fight  despite  potential  damage  suffered  in  battle.  As 
discussed  in  section  II. A.,  equipment  vital  to  a  ships  safety 
and  combat  effectiveness  is  supplied  both  normal  and  alternate 
sources  of  power.  However,  for  some  sensitive  equipment,  even 
the  brief  time  required  to  shift  from  normal  to  alternate 
power  sources  may  be  damaging.  Local  sources  such  as 
uninterruptable  power  supplies  ,  UPS's,  are  needed  during  the 
time  required  to  switch  to  an  alternate  power  source.  This 
capability  is  known  as  fight  through. 

Jb.  Reconfiguration  Techniques 

Often  there  are  several  routes  through  which  power 
may  be  supplied  to  a  vital  piece  of  equipment.  A  system  that 
automatically  reconfigures  the  distribution  system  following 
a  casualty  or  damage  might  save  a  ship  in  battle.  This  is  a 
system  level  consideration  as  compared  to  the  component  level 
solution  based  on  ABT's  and  MBT's. 
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c.  Power  Conditioning 


As  stated  in  the  Chapter  I  section  A,  advances  in 
power  electronics  are  a  driving  factor  in  the  development  of 
advanced  power  systems .  The  development  of  high  power  devices 
and  advances  in  the  control  of  solid  state  converters  has 
resulted  in  a  larger  proportion  of  the  total  electrical  load 
being  dedicated  to  solid  state  converters.  This  in  turn  has 
increased  problems  due  to  the  harmonic  frequencies  generated 
by  these  devices.  This  problem  will  most  likely  grow  worse 
over  time.  As  will  be  shown,  a  cornerstone  of  the  I  ED  system 
is  the  use  of  very  high  power  solid  state  switching 
converters . 

Partial  relief  from  this  problem  may  result  from 
improved  algorithms  for  controlling  these  converters.  It  is 
likely,  however,  that  active  power  conditioning  will  also  be 
required  to  ensure  proper  operation  of  equipment  sensitive  to 
harmonic  disturbances. 

d.  Embedded  Intelligence 

The  use  of  intelligent  control  systems  using  state 
of  the  art  monitoring  and  data  acquisition  may  allow  active 
load  control  and  configuration-based  protection  to  be 
implemented.  Active  load  control  will  ensure  that  a  system  is 
not  overloaded  by  having  too  many  devices  starting  during  a 
given  time  period  (for  example  following  a  casualty  or  when 
configuring  for  battle)  or  during  rapid  speed  changes  which 
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might  stress  an  IED  system.  The  current  solution  to  these 
problems  is  over  design  of  the  system  to  account  for  startup 
transients  and  the  use  of  LVP  controllers  as  discussed 
previously  to  ensure  time  separation  in  restarting  equipment 
following  a  casualty. 

e.  Zonal  Distribution 

Figure  3  shows  the  simplified  oneline  diagrams  of 
the  60  Hz  distribution  systems  for  a  Navy  guided  missile 
destroyer  [3]  .  The  top  diagram  shows  the  current  design  which 
matches  Fig.  2.  The  bottom  diagram  shows  the  proposed  zonal 
distribution  system.  Among  the  benefits  of  this  new  proposed 
system  are  reduced  bulkhead  penetrations  resulting  in  greater 
watertight  integrity  between  zones  and  elimination  of  over 
19000  feet  of  cable  resulting  in  a  weight  reduction  of  over  31 
tons  [3]  . 

2 .  Propulsion  Derived  Ship 1 s  Service 

By  definition  an  IED  ship  will  require  propulsion 
derived  ship's  service  (PDSS) .  PDSS  describes  a  system  in 
which  ships  service  electrical  power  is  generated,  at  least  in 
part,  from  a  ships  main  engines.  Energy  may  be  transferred 
from  the  propulsion  system  either  through  electrical  or 
mechanical  means.  Some  version  of  PDSS  with  mechanical  drive 
will  most  likely  be  an  interim  step  towards  a  complete  IED 
system.  Figures  4  and  5  show  four  potential  methods  of 
implementing  PDSS  [6] . 
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Each  of  the  four  methods  shown  requires  a  power 
converter  prior  to  the  ships  service  bus.  In  a  conventional 
gas  turbine  power  system  the  gas  turbine  generator  (GTG)  is 
speed  regulated  to  maintain  constant  frequency.  Propulsion 
turbine  speed  is  not  constant  therefore  a  variable  speed 
constant  frequency  (VSCF)  system  must  be  implemented.  Figure 
6  is  a  block  diagram  of  a  PDSS  VSCF  generating  system  [6] . 
Figure  7  is  a  100  KW  cycloconverter  tested  as  part  of  a  VSCF 
hardware  model  [6] .  A  similar  topology  might  be  appropriate 
for  shipboard  use. 

There  are  several  benefits  inherent  to  PDSS  including 

[8]  : 

•  Fuel  savings  due  to  incremental  efficiency  of  propulsion 
turbines  over  smaller  turbine  generator  sets 

•  Increased  mission  range  due  to  decreased  fuel  consumption 

•  Increased  space  due  to  removal  of  turbine  generators  and 
associated  auxiliary  equipment 

•  Arrangement  flexibility  due  to  a  more  compact  engineering 
plant 

•  Enhanced  survivability  due  to  removal  of  turbine  generator 
sets 

•  Reduced  maintenance  due  to  less  turbine  equipment 

C.  MODELLING  REQUIREMENTS  NECESSARY  TO  SUPPORT  ANTICIPATED 
DEVELOPMENTS 

As  previously  stated,  advances  in  the  design  of  ship's 
power  systems  require  extensive  testing  and  simulation  to 
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Figure  6  Block  diagram  for  propulsion  derived  ship's 
service  using  variable  speed  constant  frequency  techniques 
rS] 
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Figure  7  100  Kw  cycloconverter  basic  power  circuit  [6] 
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characterize  and  validate  system  and  component  performance. 
These  simulations  must  be  able  to  accurately  depict  system 
performance  under  all  operating  conditions  including  normal 
and  casualty  situations.  Some  examples  of  the  operating 
conditions  to  be  modelled  are: 

•  Balanced  and  unbalanced  electrical  faults 

•  Dynamic  breaking  in  an  I ED  implementation 

•  Overspeed  conditions 

•  Protective  device  trips 

•  Abrupt  speed  changes  or  shaft  reversal  (crash  back) 

Two  levels  of  detail  are  of  interest:  full  order  (or 
detailed)  and  reduced  order  simulations.  Applications  of  full 
order  simulations  include: 

•  waveform  level  studies 

•  support  of  component  development 

•  power  continuity  studies 

•  verification  of  reduced  order  models 

These  simulations  must  be  fast,  yet  they  must  retain 
electrical  transient  behavior.  Simulations  requiring  less 
detail  may  use  reduced  order  modelling  which  provides  the 
envelope  of  the  system  performance.  This  approach  is  suited 
to  simulating  electro-mechanical  transients  or  to  simulating 
electrical  systems  where  waveform  level  detail  is  not 
required.  Applications  of  reduced  order  modelling  include: 
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•  dynamic  stability 

•  control  design 

•  balanced  fault  studies 

D.  CURRENT  RESEARCH  IN  SUPPORT  OF  ANTICIPATED  DEVELOPMENTS 

1.  Machinery  Simulation  Laboratory  at  David  Taylor 
Research  Center 

The  machinery  simulation  laboratory  (MSL)  at  DTRC  is 
intended  to  provide  a  scale  model  to  test  new  power  technology 
for  shipboard  use.  While  we  will  introduce  the  MSL  at  this 
time  it  is  described  in  detail  in  the  following  chapter.  The 
laboratory  consists  of  a  turbine  emulator  driving  a  3  phase 
synchronous  generator.  The  turbine  emulator  is  composed  of  a 
DC  motor  with  solid  state  armature  supply  that  may  be  operated 
in  either  current  regulated  or  speed  regulated  modes.  A 
compliant  shaft  coupling  is  used  between  the  DC  motor  and  the 
synchronous  generator  which  allows  torque  to  be  monitored. 
Transmitted  torque  is  then  compared  to  the  simulated  response 
of  the  desired  prime  mover.  Among  the  uses  of  the  MSL  is  the 
validation  of  computer  models  developed  for  the  simulation  of 
shipboard  power  systems.  [8] 

2.  Hybrid  Analog/Digital  Computer  Simulations  at  Purdue 
University 

Purdue  University  has  designed  and  extensively  used  a 
computer  based  on  a  combination  of  analog  and  digital 
technology  [8]  ,  [9]  .  This  system  has  been  used  to  conduct 
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system  studies  for  the  Navy.  Among  the  benefits  of  this 
approach  is  speed.  The  speed  of  the  simulation  is  independent 
of  the  system  complexity  due  to  the  use  of  parallel 
integration.  Purely  digital  techniques  use  serial  integration 
which  creates  a  direct  link  between  solution  speed  and  system 
complexity.  The  drawbacks  of  an  analog  simulation  include 
high  overhead  in  time  and  manpower  setting  up  a  simulation  and 
a  limitation  on  the  size  of  the  system  to  be  modelled  based  on 
equipment  available.  [6] 

3.  Digital  Computer  Modelling 
a.  Purdue  University 

Other  work  at  Purdue  University  consists  of  using 
the  commercial  package  Advanced  Computer  Simulation  Language 
(ACSL)  to  model  systems,  subsystems,  and  components  related  to 
ships  power  systems.  This  effort  uses  strictly  digital 
simulation  techniques. 

Further  work  at  Purdue  is  in  the  area  of  reduced 
order  modelling.  Ship's  power  systems  are  "stiff",  that  is 
various  time  constants  within  the  system  are  widely  separated. 
The  small  time  step  required  for  the  numerical  solution  of  the 
fast  transients  makes  it  infeasible  to  conduct  the  long 
simulations  required  to  model  the  mechanical  transients 
associated  with  an  I ED  system.  If  one  can  neglect  the  fast 
transients  then  the  solution  becomes  much  faster.  [6] 
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b.  Massachusetts  Institute  of  Technology 

The  Massachusetts  Institute  of  Technology  has 
addressed  the  problems  associated  with  modelling  ships  power 
systems  with  an  ambitious  project  to  develop  a  new  software 
package,  WAVESIM.  WAVESIM  addresses  the  stiffness  problem 
associated  with  ships  power  systems  through  the  use  of 
waveform  representation  of  system  variables.  The  use  of 
waveform  techniques  allows  much  larger  time  steps  than  do 
traditional  integration  schemes.  A  detailed  explanation  of 
the  methods  employed  by  WAVESIM  is  given  in  Chapter 
V.  [1]  [4]  [12]  [13] 
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III.  REQUIRED  SOFTWARE  CAPABILITIES 


An  evaluation  requires  a  set  of  metrics.  This  paper 
proposes  a  set  of  metrics  to  be  used  for  the  evaluation  of 
detailed  waveform  analysis  tools  used  to  simulate  shipboard 
power  systems.  General  requirements  for  a  tool  of  this  nature 
are  given  by  [1]  ,  [4]  ,  and  [13]  .  The  following  paragraphs 
combine  and  extend  the  requirements  set  forth  in  [1]  ,  [4] ,  and 
[12]  . 

A.  SYSTEM  AND  COMPONENT  LEVEL  MODELLING  ENVIRONMENT 

The  proposed,  system  simulator  must  be  modular,  allowing 
for  the  interconnection  of  various  component  models  into 
widely  varying  system  topologies.  The  ability  to  model  highly 
non-linear  components  and  components  that  introduce 
discontinuities  to  the  system  is  essential . 

Available  component  models  should  include  at  a  minimum: 
prime  movers,  synchronous,  induction,  and  dc  machines;  solid 
state  converters;  and  various  passive  load  representations. 

The  simulation  should  be  intuitive  in  that  both  component 
models  and  system  models  should  reflect  their  physical 
counterparts . 

Lastly,  it  is  essential  that  the  user  be  able  to  add 
component  models.  Component  model  addition  should  involve  a 
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minimum  of  effort  in  the  transition  from  physics  to  computer 
code. 

B.  ROBUSTNESS  WHEN  SIMULATING  NON-LINEAR  AND  RAPID  SWITCHING 
TOPOLOGIES 

The  simulation  must  be  able  converge  to  the  solution 
representing  the  actual  operation  of  the  component  devices. 
A  simple  example  is  to  disallow  currents  from  existing  in  the 
reverse  bias  direction  of  a  diode.  Models  must  be  capable  of 
detailed  representation  of  transient  behavior  during  and 
following  topology  changes  due  to  events  such  as  solid  state 
switching,  circuit  breaker  closing,  faults,  or  equipment 
start -up/ shut  down. 

C.  CORRECTNESS  OF  SOLUTION 

The  simulation  must  accurately  depict  the  characteristics 
of  interest  in  support  of  detailed  waveform  analysis.  Sources 
of  numerical  error  must  be  addressed,  quantified  and  bounded. 
Some  form  of  feedback  as  to  the  reliability  of  the  solution 
must  be  provided  by  the  simulation.  Adequate  validation  must 
be  performed  against  real  world  data  or  accepted  similar 
simulations. 

D.  SOFTWARE  DOMAIN 

A  means  of  importing  field  test  data  should  be  provided. 
The  ability  to  plot  imported  data  superimposed  on  simulation 
results  provides  a  valuable  tool.  Further,  a  means  of  using 
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imported  field  test  data  to  drive  component  models  is 
desirable.  Models  should  be  developed  such  that  various 
control  parameters  are  easily  changed  to  aid  in  control  system 
synthesis . 

E.  IMPLEMENTATION  OPTIONS 

The  simulation  must  be  fully  portable.  The  chosen 
programming  language,  algorithms,  and  documentation  methods 
should  be  written  to  provide  long  term  maintainability. 
Graphics  output  presents  particular  problems  with  portability. 
One  way  to  preclude  problems  in  this  area  is  to  ensure 
graphics  output  is  handled  by  a  separate  and  isolated 
subprogram.  Data  input  and  output  options  should  include  both 
interactive  and  file  modes. 

F.  EASE  OF  USE 

The  simulation  should  be  user  friendly.  Both  interactive 
and  file  input  modes  should  be  intuitive.  Display  of 
simulation  results  should  be  flexible.  The  capability  to 
easily  plot  potentials  with  respect  to  a  reference  or  across 
devices,  to  plot  flow  variables  and  to  manipulate  data 
(shifting,  multiplication,  addition,  subtraction,  integration, 
and  differentiation)  are  required.  Simulation  results,  in  the 
form  of  time  series  data,  should  be  saved  to  a  data  file  in 
standard  (ASCII)  format.  The  graphics  subprogram  should  be 
able  to  present  previously  saved  data  from  these  files. 


25 


6.  SOFTWARE  SPEED  VERSUS  SYSTEM  COMPLEXITY 


Due  to  the  potential  complexity  of  the  systems  to  be 
modeled,  the  simulation  must  be  fast  while  retaining  transient 
behavior.  For  specialized  studies,  reduced  order  models 
should  be  implemented  to  enhance  simulation  speed  while 
maintaining  accurate  representation  of  the  characteristics  of 
interest . 

H.  CONTINUED  SUPPORT 

A  program  applied  to  ship  power  system  design  and  analysis 
will  be  used  in  arriving  at  decisions  involving  millions  of 
dollars  in  government  expenditures  and  may  be  used  over  a 
period  of  many  years.  It  is  essential,  therefore,  that  such 
a  program  carry  with  it  a  long  term  commitment  of  support. 
Essentially  this  requires  a  package  to  be  commercially 
supported  or  supported  by  an  organization  equipped  to  make 
this  kind  of  commitment. 
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IV.  DETAILED  MODELS  NECESSARY  TO  SIMULATE  SHIPS  POWER  SYSTEMS 


The  required  models  are  placed  into  two  groups.  Those 
models  required  to  simulate  the  MSL  are  given  first.  These 
models  will  be  used  to  demonstrate  WAVESIM.  The  second  group 
of  models  are  derived  from  studies  conducted  by  P.C.  Krause 
for  use  in  work  conducted  for  the  Navy  by  Purdue  University 
and  PC  Krause  and  Associates  [9] .  While  the  models  included 
are  not  all-inclusive  they  provide  a  representative  sample  of 
the  type  of  models  that  are  required.  The  only  models  covered 
in  detail  are  those  used  in  the  testing  of  WAVESIM. 

A.  MACHINERY  SIMULATION  LAB  MODELS 

A  block  diagram  of  the  machinery  simulation  lab  at  David 
Taylor  Research  Center  is  given  in  Fig.  8.  This  lab  consists 
of  a  DC  motor  which  operates  in  either  a  speed  regulated  or 
current  regulated  mode  driving  a  synchronous  3  phase 
generator.  The  DC  motor,  current  regulator,  and  speed 
regulator  make  up  the  turbine  emulator  which  provides  dynamic 
emulation  of  various  gas  turbine  and  diesel  prime  movers.  [8] 
1 .  Speed  Regulator 

The  block  diagram  for  the  speed  regulator  is  given  in 
Fig.  9.  In  this  model  a  reference  speed,  wref,  is  compared  to 
the  actual  motor  speed,  wm,  resulting  in  the  output  of  a 
reference  current,  iref,  which  goes  to  the  current  regulator. 
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Figure  8  Block  Diagram  of  the  machinery  simulation  lab  at 
David  Taylor  Research  Center  [8] 
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Figure  9  Block  diagram  for  machinery  simulation  lab  speed 
regulator  [8] 
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When  used  to  emulate  a  prime  mover,  the  reference  speed  is 
generated  by  a  comparison  of  the  transmitted  shaft  torque 
between  the  DC  motor  and  synchronous  generator  to  simulate  the 
response  of  the  prime  mover  of  interest. 

The  variable  LSRA  in  Fig.  9  represents  the  amplifier 
limits  using  a  non  wind-up  limiter  as  described  in  [13] .  This 
limiter  is  modelled  simply  as  a  minimum/maximum  filter  passing 
the  actual  output  value  of  a  block  if  the  output  is  within  the 
limits  and  passing  either  the  minimum  or  maximum  value  of  the 
filter  if  the  output  of  the  block  is  outside  these  limits. 

2 .  Current  Regulator 

The  current  regulator  controls  the  motor  current  by 
comparing  the  reference  current  provided  by  the  speed 
regulator,  iref,  to  the  motor  current,  im.  Figure  10  is  a 
block  diagram  of  the  current  regulator.  Note  the  use  of  the 
non- windup  limiter  once  again.  The  output  of  the  current 
regulator  is  the  DC  motor  armature  voltage,  vm. 

3 .  DC  Motor  Model 

Figure  11  shows  a  simple  first  order  model  chosen  to 
model  the  DC  motor.  A  more  detailed  model  would  include  core 
losses,  windage  and  friction  losses  as  detailed  in  [14] .  The 
error  generated  by  this  approximation  as  reported  in  [8]  is  on 
the  order  of  3.5  percent  at  rated  conditions. 
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Figure  11 

simulation 


Block  diagram  of  DC  motor  of  the  machinery 
lab  [8] 
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4.  Three  Phase  Synchronous  Generator  Model 

The  synchronous  generator  model  used  to  illustrate  the 
capabilities  of  WAVESIM  is  derived  in  [15]  .  This  model 
consists  of  seven  first  order  differential  equations  given  in 
the  rotor  reference  frame. 

The  equations  describing  the  electrical  and  mechanical 
behavior  of  the  three  phase  synchronous  generator  are: 

-5-f  =  -bt*  +  y2-  +  +  «0VSin«  (1) 

ac  1ad  Aad 


=  “  IjT3-  ~  IS—  *  -m— 


+  WgVCOSfi 
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The  following  is  a  partial  list  of  the  constants  and  states 
found  in  (1)  to  (7)  : 

•  7d  Per-unit  direct  axis  flux  linkage 

•  7q  Per-unit  quadrature  axis  flux  linkage 

•  eg"  Quadrature -axis  voltage  behind  subtransient  reactance 

•  ed"  Direct -axis  voltage  behind  subtransient  reactance 

•  eqn  Quadrature -axis  voltage  behind  transient  reactance 

•  <i >0  Base  frequency 

•  6  Rotor  phase  angle  with  respect  to  synchronous 

reference 

•  do "  Direct -axis  open- circuit  subtransient  time  constant 

•  Tqo"  Quadrature -axis  open-circuit  subtransient  time 

constant 

•  Tdo'  Direct-axis  open-circuit  transient  time  constant 

•  H  Ratio  of  mechanical  energy  at  rated  speed  to  base 

power 

•  Xq  Quadrature -axis  synchronous  reactance 

•  Xq”  Quadrature -axis  subtransient  reactance 

•  Xd'  Direct-axis  transient  reactance 

•  Xd”  Direct-axis  subtransient  reactance 

Depending  on  the  nature  of  the  simulation  to  be  conducted, 
either  the  phase  currents  or  phase  voltages  are  assumed  known. 
For  balanced  three  phase  operation,  the  phase  quantities  may 
be  mapped  to  or  from  the  rotor  reference  frame  using  Parks 
transformation: 
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A  detailed  explanation  of  these  equations  is  beyond 
the  scope  of  this  thesis.  An  excellent  short  introduction  to 
modelling  synchronous  machines  is  given  in  [15] .  Analysis  of 
electric  machinery  with  thorough  treatment  of  reference  frame 
theory  and  Park's  transformation  is  presented  in  [16]  and 
[17]  . 

5.  Voltage  Regulator  Model 

Figure  12  shows  the  block  diagram  of  the  excitation 
system  for  the  MSL.  The  voltage  regulator  is  of  Type  II  as 
defined  by  the  IEEE  [18]  .  Vref  is  a  reference  voltage  supplied 
to  the  regulator  while  Vt  is  the  average  terminal  voltage.  The 
second  half  of  the  block  diagram  accounts  for  non-linear 
saturation  effects.  The  output  of  this  element  is  the  field 
excitation  voltage,  exfd. 
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Figure  12  Block  diagram  of  the  voltage  regulator  for  the 
machinery  simulation  lab  [9] 
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6.  Resistive -Inductive  Load 

For  the  purpose  of  testing  WAVESIM  only  a  wye 
connected  R-L  load  model  is  available.  The  phase  voltages  are 
inputs  and  the  phase  currents  are  outputs.  The  following 
equations  represent  a  wye  connected  RL  load  [8] : 


VL  = 


VL  = 


VG  J-Cl 


diaL 

“b 

dt 

* 

dibL 

**> 

dt 

dt 

(11) 


(12) 


(13) 


B.  ADDITIONAL  MODELS  FOR  SHIPBOARD  ELECTRICAL  DISTRIBUTION 

SYSTEM 

A  specialized  shipboard  electrical  distribution  system 
studied  by  Purdue  University  is  shown  in  Fig.  13.  The  system 
consists  of  two  identical  ships  service  turbine  generators 
(SSTGs)  supplying  main  and  vital  loads.  Vital  loads  are 
supplied  alternate  power  from  a  battery  backup.  DC/AC  and 
AC/DC  conversion  is  accomplished  via  ships  service  motor 
generator  (SSMG)  sets  consisting  of  a  synchronous  AC  machine 
mechanically  coupled  to  a  DC  machine.  Normal  operation  is  for 
the  AC  machine  to  act  as  a  motor  driving  the  DC  machine  to 
charge  the  batteries.  In  a  casualty  the  batteries  would 
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supply  the  DC  machine  acting  as  a  motor  driving  the  AC  machine 
to  provide  vital  60  Hz  power.  [9] 

This  system  will  not  be  simulated  in  this  thesis  since 
adequate  models  for  WAVESIM  are  not  available;  rather,  it  is 
an  illustration  of  the  type  of  study  that  a  tool  such  as 
WAVESIM  is  required  to  perform.  Some  models  required  by  this 
study  in  addition  to  those  presented  in  the  previous  section 
are: 

•  Induction  motors 

•  Solid-state  power  converters 

•  Prime  movers 

•  Static  excitation  system 

•  DC  conversion  system 


1.  Induction  Motor  Model 

"The  induction  motor  is  the  workhorse  of  the  electric 
power  industry"  [16:pl59] .  This  is  true  of  the  Navy  as  well. 
The  Navy  relies  heavily  on  induction  motors  to  drive  the 
hundreds  of  pumps,  fans,  and  compressors  found  aboard  ship. 
The  equations  governing  an  induction  machine  in  an  arbitrary 
reference  frame  and  using  per  unit  quantities  are  given  in 
[8]  . 


V  =  r  i  + 
vqs  '4-s-tqs 


+  — 
dt 


(14) 


39 


(15) 


„  r-  V  _  .  d  ♦dS 

ds  ~  rs*ds  —  *qs  -g^  — 


Vbs  -^s^Os  + 


_d  *3s 
dt 


n_  <  ,  (  ^  ~  \b  ,  d  *frqr 

„  _  ,  ,■  _f“-  “H  L  ^  d  *dr 

0  -  r^*  {— ^~J*"  Silt 


r\  -  r  j  .  &  'fcl* 

0  "  r'lar  +  dtlt 


The  per  unit  flux  linkages  are  given  by 

—  -^ls-^qs+-^  ^os+-^qr^ 


*^ds  =  -^Is-^ds  +  -^M^ds  +  -^dr^ 


'fbs 


Os 


^qr  -^Ir-^qr  +  ^qs  +  -^qr^ 
^dr  =  ■‘^Ir-^dr  +  •^M^ds+-^dr^ 
♦br  =  -^Ir  -*0r 

Finally,  the  per  unit  torque  is  given  by 


*^qr  -*dr  “  *Klr  ^ 


qr 


(16) 


(17) 


(18) 


(19) 


(20) 

(21) 

(22) 

(23) 

(24) 

(25) 


(26) 


Rather  than  detail  every  variable  name,  the  following  general 
guidelines  are  provided: 


40 


•  a  is  a  rotational  speed  in  radians  per  second 

•  f  is  the  per-unit  flux  linkage 

•  X  is  a  reactance 

•  subscript  b  base  quantity 

•  subscript  q  refers  to  quadrature  axis 

•  subscript  d  refers  to  direct  axis 

•  subscript  0  refers  to  zero  sequence  axis 

•  subscript  s  refers  to  stator  variables 

•  subscript  r  refers  to  rotor  variables 

•  subscript  1  refers  to  a  leakage  reactance 

•  subscript  M  refers  to  the  magnetizing  reactance 

In  the  above  equations  u  is  the  rotational  speed  of  the 
reference  frame.  For  a  stationary  reference  frame  a  is  set  to 
zero.  Translation  to  and  from  the  rotor  reference  frame  may  be 
accomplished  by  Park's  transformation's,  as  given  in  section 
A. 4  of  this  chapter. 

2 .  Prime  Movers 

Prime  movers  for  naval  applications  vary  from  internal 
combustion  engines  such  as  large  diesel  engines  and  gas 
turbine  engines  to  high  pressure  steam  turbines.  A  simple 
second  order  model  of  a  steam  turbine  prime  mover  which 
accounts  for  both  plant  and  servo/steam  valve  time  constants 
is  given  in  Fig.  14. 
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Figure  14  Second  order  model  of  a  steam  turbine  prime 
mover  [10] 


3.  Static  Excitation  System 

Figure  15  which  is  a  representative  static  excitation 
system  [9]  .  This  diagram  is  provided  as  an  example  of  the 
complexity  of  the  models  that  need  to  be  modelled.  A  detailed 
explanation  for  this  excitation  system  may  be  found  in  [9]  or 
[13]  . 

4 .  DC  System 

The  distribution  system  shown  in  Fig.  13  is  highly 
specialized  yet  it  demonstrates  an  important  requirement  for 
power  system  simulations:  the  ability  to  simulate  abruptly 
changing  topologies.  Figure  16  shows  the  DC  subsystem  of  Fig. 
13 .  Note  that  a  switch  determines  the  mode  of  operations  of 
the  DC  subsystem.  In  normal  operation  the  DC  machine  acts  as 
a  generator.  The  feedback  is  based  on  the  field  excitation  to 
maintain  proper  DC  voltage  output.  In  emergency  operation, 
feedback  is  based  on  a  reference  speed  necessary  to  maintain 
the  appropriate  frequency  of  the  synchronous  machine  now 
acting  as  a  generator.  [10] 

5.  Solid-state  Power  Converters 

The  three  phase  uncontrolled  diode  rectifier  shown  in 
Fig.  17  is  a  simple  three  phase  power  converter.  This  device 
produces  an  uncontrolled  DC  voltage  output  from  a  three  phase 
sinusoidal  input.  Figure  18  is  the  equivalent  circuit  for 
Fig.  17  as  given  by  [19:p.27]. 
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Figure  15  Representative  static  excitation  system  [10] 
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Figure  16  DC  subsystem  of 
backup  [10] 


a  ships  power  system  with  DC 
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V.  DETAILED  DESCRIPTION  OF  WAVESIM 


A.  WHAT  IS  WAVESIM 

WAVESIM  is  a  computer  program  developed  by  LT  Norbert 
Doerry,  USN  at  the  Massachusetts  Institute  of  Technology  [1]  . 
It  is  designed  to  perform  detailed  digital  simulation  of 
systems  of  linear  and  nonlinear  devices.  The  ultimate 
objective  in  developing  WAVESIM  is  to  simulate  ship's  IED 
systems.  With  properly  developed  device  modules,  however,  it 
should  be  capable  of  modelling  any  system  of  non-linear 
devices . 

In  its  current  form  WAVESIM  is  a  demonstration  project 
only.  It  consists  of  a  collection  of  source  code  files 
written  in  the  C  programming  language,  MATLAB  script  files, 
and  specially  formatted  files  defining  both  device 
input/output  and  system  connections.  The  types  of  files  are: 

•  WAVESIM  source  code  (.c  files  and  .h  files) 

•  Waveform  operators  (MATLAB  script  files) 

•  Definition  files  (special  format  files  with  extension 

.  def) 

•  Device  constitutive  files  (MATLAB  script  files) 

•  Input  files  (special  format  files  with  extension  .in) 
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1.  WAVESIM  Source  Code 

WAVESIM  itself  is  written  in  ANSI  standard  C  and  is 
highly  portable.  The  program  has  been  compiled  and  executed 
on  a  variety  of  platforms  including  Sun  SPARC2  work  stations, 
VAXS TAT ION  3100  and  VAX  11/785,  and  IBM  Personal  Computers. 
The  output  of  WAVESIM  is  a  simulation  file  in  the  form  of  a 
script  or  .m  file  for  execution  under  MATLAB  by  Mathworks, 
Inc.  WAVESIM  takes  advantage  of  MATLAB' s  extensive  vector 
mathematics  capability.  The  simulation  files  are  intended  to 
be  used  as  provided.  The  user  should  have  no  reason  to  modify 
tnem. 

2 .  Waveform  Operators 

WAVESIM  uses  waveform  representation  of  interface 
variables  and  gives  the  user  the  option  of  several 
representations.  The  waveform  operators  supplied  with  the 
program  as  MATLAB  functions  perform  all  the  required 
mathematical  operations.  These  operations  include  converting 
from  one  representation  to  another,  finding  derivatives, 
computing  integrals  over  an  interval,  smoothing,  and 
arithmetic  operations.  These  waveform  operators  are  used  both 
by  the  device  constitutive  files  and  by  the  simulation  files 
generated  by  WAVESIM.  These  files  are  intended  to  be  used  as 
provided.  The  user  should  have  no  reason  to  modify  them. 
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3.  Definition  Files 

Files  describing  the  device  have  the  extension  .def. 
Each  device  must  be  defined  in  the  special  format  of  the  .def 
file.  The  device  definition  file  provides  the  following 
information  [1] : 

•  Name  of  device  type 

•  Number,  names,  and  default  values  of  parameters 

•  Number,  names,  and  default  initial  values  of  states 

•  Number  of  terminals 

•  Terminal  names 

•  Terminal  Types  (normal  or  information) 

•  Flow  variable  types  (import  or  export) 

•  Potential  variable  types  (import  or  export) 

•  Terminal  KCL  group  numbers 

•  Device  structural  jacobian 

Users  desiring  to  add  devices  to  the  library  must 
prepare  a  new  name. def  file  in  the  format  given  in  [l] .  The 
main  .def  file  is  device. def.  This  file  contains  debugging 
information  and  lists  all  defined  devices  either  with  complete 
definitions  or  with  include  statements.  Device. def  must  be 
present  for  WAVESIM  to  operate.  New  device  definitions  may  be 
added  directly  to  device. def  or,  the  preferred  method  is  to 
develop  a  separate  name. def  file  and  add  it  to  the  library 
using  an  include  statement  in  device. def. 


50 


The  information  provided  by  the  device  definition  file 
is  sufficient  to  combine  devices  and  reduce  the  resulting 
system  of  equations  into  solvable  blocks.  However,  no 
information  is  provided  on  the  constitutive  relationship  of 
the  systems  variables.  This  information  is  provided  by  the 
device  .m  files. 

4.  Device  Files 

Device  files  are  MATLAB  script  files  that  detail  the 
constitutive  relationships  of  the  import  and  export  variables 
of  a  device.  Doerry  lists  the  properties  of  these  files  in 
terms  of  resources  and  products  of  the  file  [1]  .  The 
resources  are: 

•  Waveform  type 

•  Import  variable  waveform 

•  Parameter  values 

•  Value  of  states  at  beginning  of  time  interval 

•  Time  Structure 

•  Beginning  time  of  interval 

•  Ending  time  of  interval 

•  Minimum  time  of  interval  of  interest 

•  Continuation  parameter 

With  these  resources  WAVESIM  uses  the  device  file  to  calculate 
the  following  products: 

•  Export  variable  waveforms 
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•  Device  jacobian  matrix 

•  Value  of  states  at  end  of  time  interval 

•  Recommended  time  structure 

•  Recommended  recalculation  time  this  interval 

•  Recommended  ending  time  of  next  interval 

5.  Input  Files 

Input  files  have  the  extension  .in  and  are  written  by 
the  user.  The  input  files  initialize  WAVESIM  in  regard  to  the 
system  being  studied.  The  file  is  command  oriented  with  seven 
available  commands  [1] : 


•  Debug 

•  Default 

•  Device 

•  Include 

•  Node 

•  Time 


Print  debug  information 
Set  default  system  parameters 
Specify  device  information 
Include  another  file 
Specify  node  parameters 
Control  time  increments 


Commands  may  be  either  a  single  line  or  they  may  be  multiple 
lines  as  designated  by  placing  the  command  on  a  line  by  itself 
followed  by  a  list  of  subordinate  commands  and  ending  with  the 
word  END. 

The  following  is  a  brief  summary  of  each  of  the  commands, 
a.  Debug 

WAVESIM  has  numerous  built  in  debugging  modules. 
These  modules  print  out  various  simulation  data  as  the 
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simulation  progresses  to  allow  the  user  to  identify  problems 
with  the  simulation.  The  type  of  information  available 
includes  all  the  devices  defined  by  WAVESIM,  structural 
jacobian,  file  data,  and  construction  of  various  blocks. 
b.  Default 

Control  of  the  simulation  is  exercised  with  the 
default  command.  In  all  there  are  eleven  subordinate  default 
commands,  most  with  multiple  elements.  The  subordinate 
commands  are: 

•  alpha- -Control  the  continuation  parameter. 

•  check- -Error  checking  flags. 

•  diverge- -Control  of  divergence  test. 

•  error- -Set  default  error  levels. 

•  gmin--Set  default  node  leakage  conductance. 

•  max- -Set  maximum  number  of  iterations. 

•  nbr- -Control  number  of  coefficients  in  waveform 

representation . 

•  rmin--Set  default  series  resistance. 

•  scale- -Set  default  node  scaling  factors. 

•  waveform  content- -Set  waveform  content  limits. 

•  wtype--Set  default  waveform  type. 

With  experience,  tight  control  of  the  simulation  and  error 
control  are  achievable  using  the  default  settings. 
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c.  Device 


The  topology  of  the  system  being  modelled  is 
specified  by  the  device  command  and  its  subordinates.  A  given 
system  may  contain  multiple  occurrences  of  a  given  device  type 
(i.e.  multiple  resistors) .  Each  particular  device  is  given  a 
unique  name  by  which  it  is  identified  throughout  the 
simulation.  The  name  is  assigned  as  part  of  the  device 
command.  The  device  subordinate  commands  are: 

•  Terminal 

•  Parameter 

•  States 


(1)  Terminal 

The  terminal  subordinate  command  assigns 
terminals  to  specific  nodes  and  hence  specifies  the  system 
topology.  Each  terminal  of  a  device  must  be  assigned  to  a 
node.  The  rules  governing  assignment  of  terminals  to  nodes 
are  given  in  [1]  . 

(2)  Parameter 

The  use  of  parameters  in  device  definitions 
allows  for  one  device  model  to  be  written  for  each  device  type 
For  example,  a  synchronous  generator  may  be  a  2  MW  unit  with 
its  associated  parameters  or  it  may  be  a  750  KW  unit  with  its 
associated  parameters  but  both  may  use  the  same  device  model. 
Default  values  of  parameters  are  specified  in  the  device. def 
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files. 


Parameters  in  the  device  command  override  default 


values.  (3)  States 

Doerry  gives  the  following  definition  of 


states  [1:  p.46] : 

STATES  are  variables  whose  values  are  stored  for  a  given 
time  for  later  use.  States  can  be  used  for  example,  to 
store  the  constant  of  integration  for  a  dynamic  equation. 
States  can  also  be  used  to  store  the  operating  mode  for  a 
given  device.  In  general,  if  the  value  of  a  given 
variable  depends  on  the  previous  value  of  another  variable 
that  other  variable  is  a  state. 

The  device  states  subordinate  command  allows  the  user  to 

initialize  a  given  state.  This  may  be  crucial  to  the 

simulation  since  an  iterative  technique  is  used  to  solve  the 

systems  or  algebraic  equations  representing  the  plant.  The 

iterative  method  used  is  Newton -Raphson.  A  region  of 

convergence  is  located  around  each  solution  to  a  non-linear 

algebraic  equation.  Proper  initializing  of  the  system  may  be 

key  to  obtaining  the  correct  solution.  This  is  discussed 

further  in  section  B.4  of  this  chapter. 

d.  Include 

The  include  command  inserts  a  specified  text  file 
at  the  location  of  the  include  statement. 

e.  Node 

The  node  command  specifies  node-specific 
parameters.  For  instance,  the  user  may  specify  a  resistance 
(Rmin)  or  conductance  (Gmin)  to  ground  to  reduce  linear 
dependencies  among  the  systems  of  equations  at  a  particular 
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node.  Doerry  details  the  use  of  Rmin  and  Gmin  to  address 
singularities  in  the  jacobian  matrix  [l:p.53] .  Maximum 
equation  errors  and  maximum  variable  corrections  may  also  be 
specified  for  a  given  node.  Other  subordinate  node  commands 
include  name,  to  name  a  particular  node,  and  scale,  to 
override  default  scale  factors  for  potential  or  flow  variables 
at  a  particular  node. 
f .  Time 

The  time  command  specifies  start  and  finish  times 
for  the  simulation  as  well  as  allowing  the  insertion  of  break 
points  and  controlling  time  increments.  Break  points  can  be 
inserted  at  the  time  of  known  discontinuities  to  speed  the 
simulation  by  ensuring  that  a  waveform  boundary  occurs  at  the 
discontinuity.  The  problem  of  waveform  representation  at 
discontinuities  and  the  difficulties  this  imposes  when 
conducting  a  detailed  simulation  involving  power  converters  is 
addressed  in  section  B.3  of  this  chapter. 

Time  increment  control  includes  the  ability  to 
specify  maximum,  minimum,  optimum,  and  initial  time  increment 
as  well  as  the  minimum  time  interval  of  interest  which  is 
used  to  by  devices  to  smooth  export  variables  and 
discontinuities.  WAVESIM  varies  the  simulation  time  for  each 
step  depending  on  the  progress  of  the  simulation. 
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g.  Plot 


The  plot  command  specifies  which  system  variables 
need  to  be  converted  to  data  series  and  plotted  in  MATLAB. 
Only  variables  specified  by  the  plot  command  are  converted  to 
data  series  for  output. 

B.  SOLUTION  METHOD  USED  IN  WAVES IM 
1.  Device  Modelling 

WAVE SIM  device  models  are  developed  using  terminal 
descriptions.  With  terminal  descriptions  a  device  is  assigned 
variables  for  flow  and  potential  at  each  terminal.  Device 
constitutive  relationships  are  then  developed  in  terms  of 
these  variables.  In  contrast,  most  engineers  are  familiar 
with  branch  descriptions  in  which  flow  variables  are  written 
in  terms  of  the  relative  potential  difference  between 
terminals.  Figure  19  is  a  simple  example  of  branch  vs 
terminal  descriptions  for  a  two  terminal  device. 

Doerry  defines  the  variables  v,,  i,,  v2,  and  i2  in  Fig. 
19  as  interface  variables.  These  variables  are  the  means  by 
which  various  devices  interact.  Interface  variables  may  be 
either  potential  variables  (voltage,  speed,  pressure,  etc.)  or 
flow  variables  (current,  torque,  fluid  flow  rate,  etc.). 
Interface  variables  are  further  characterized  as  imports, 
which  devices  see  as  a  resource,  or  as  exports,  which  the 
device  sees  as  a  product.  [I:p44] 
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Branch  Terminal 

description  description 

Figure  19  Branch  vs  Terminal  Description  [l:p.l5] 

Doerxy  describes  the  device  description  in  two  ways 

...  an  organized  manner  for  describing  the  characteristics 
of  a  physical  component  [l:p.44] ; 

...  a  means  for  generating  export  variables  based  on  the 
values  of  the  import  variables,  states,  parameters, 
continuation  parameter,  and  time  [l:p.45]. 

The  continuation  parameter  associated  with  the  Newton- Raphson 

method  of  solving  systems  of  non-linear  algebraic  equations 

will  be  discussed  in  part  4  of  this  section. 

For  the  device  shown  in  Fig.  19  assume  the  flow 

variables,  and  i2,  are  defined  as  export  variables  and  the 

potential  variables,  V!  and  v2,  are  defined  as  import 
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variables.  The  device  constitutive  equations  might  then  be 
expressed  as 


v2) 

(27) 

4=^2  irV-2) 

(28) 

These  equations,  expressed  on  the  device  level,  are  used  by 

WAVESIM  to  satisfy  the  system  equations. 

Two  more  products  of  the  device  object  are  the  device 

jacobian  matrix  and  the  device  structural  jacobian  matrix. 

The  device  jacobian  matrix  expresses  the  sensitivity  of  a 

given  export  variable  to  a  change  in  a  given  export  variable. 

The  device  structural  jacobian  matrix 

.  .  .  describes  the  properties  of  the  elements  of  the  device 
jacobian  matrix  for  a  given  type  of  variable 
representation  without  providing  any  values  [l:p.48] 

For  an  ordered  set  of  import  variables 


(29) 
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and  an  ordered  set  of  export  variables 


the  device  jacobian  matrix  is  given  by 
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(31) 


The  elements  of  the  device  jacobian  matrix  may  be  scalar 
values  or  matrices  depending  on  how  variables  are  represented 
(i.e.  by  real  numbers  or  by  waveform  representation) . 

[1:  pp. 47-48] 

The  device  structural  jacobian  matrix  gives  the 
mathematical  nature  of  the  device  jacobian  matrix.  Table  II 
lists  the  codes  used  in  the  device  structural  jacobian. 
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TABLE  II  CODES  FOR  STRUCTURAL  JACOBIAN  MATRIX  [1:  p.  48] 

Code 

Type  of  Matrix 

0 

Zero  matrix  (all  elements  always  zero) 

I 

Identity  matrix  (always  the  identity  matrix) 

D 

Diagonal  (always  a  linear  main  diagonal  matrix) 

L 

Linear  matrix  (  The  elements  are  always  constant . ) 

A 

Nonlinear  AC  matrix  (See  note  below.) 

N 

Nonlinear  Matrix  (The  elements  may  not  be  constant.) 

U 

Unknown  (The  dependence  is  unknown.  Treat  as  nonlinear.) 

Doerry  defines  a  nonlinear  AC  matrix  as 

.  .  .  one  for  which  the  constant  component  of  the  export 
variable  depends  only  on  the  constant  component  of  the 
import  variable.  The  other  components  of  the  export 
variable  cannot  depend  on  the  constant  component  of  the 
import  variable  but  are  not  restricted  in  any  other  way 
[l:p.  48]  . 


Device  structural  jacobians  are  combined  by  WAVESIM  to  build 
a  system  structural  jacobian  which  will  be  discussed  in  the 
following  section.  Further,  the  device  structural  jacobian 
identifies  which  elements  must  be  recalculated  between 
iterations  if  an  iterative  type  solution  is  used. 
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Lastly,  device  models  provide  WAVESIM  with  a 
recommended  recalculation  time  for  the  current  interval. 
Waveform  representation  of  variables  may  be  inaccurate  at 
discontinuities  unless  a  waveform  boundary  falls  at  the  time 
of  the  discontinuity.  Therefore,  if  a  device  constitutive 
equation  is  discontinuous,  the  model  recommends  a  new  interval 
causing  the  waveform  boundary  to  fall  at  the  discontinuity. 
WAVESIM  will  generally  use  the  minimum  recalculation  time 
suggested  by  any  of  the  devices . 

2 .  Building  a  System  of  Equations  and  Block  Reduction 

WAVESIM  uses  an  adaptation  of  modified  nodal  analysis 

to  develop  system  equations.  Figure  20  is  used  by  Doerry  to 

demonstrate  this  procedure. 

System  variables  are  defined  as 

The  minimum  set  of  variables  from  which  all  of  the  device 
import  and  export  variables  can  be  derived  [l:p.51] . 

System  variables  that  must  be  solved  for  in  Fig.  20  are  is1, 

ici#  ©o#  ©i»  and  e2.  This  information  is  based  on  a  knowledge 

of  how  the  particular  devices  in  Fig.  20  are  modelled:  which 

terminal  variables  each  model  uses  as  imports  and  which  are 

exports.  The  KCL  equations  are  written  for  each  node  using 

the  constitutive  equations  yielding: 


^i+STrjri  (e,,e2)  =0 

(32) 

■^ci+9rR.iR2  (e, ,  e2)  =0 

(33) 
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Figure  20  Terminal  description  of  an  RC  circuit  [l:p.23] 

iG+9rs.,s2(isi»eo)  +9cjc2(ici'eo)  =0  (34> 

Notice  that  the  device  constitutive  equations  are  not  known  on 
the  system  level.  Next  the  potential  difference  equations  are 
written  between  each  node  potential  and  each  export  potential 
variable  associated  with  a  normal  or  information  node. 

ei_STc.vci  (ici '  eo)  =0  (35) 

e2-grs_vsi  <^si ,  e0)  =0  (36) 

eo-sr(LvG  <  *£)  =0  (37> 
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The  total  number  or  system  equations  to  be  solved  is  equal  to 
twice  the  number  of  normal  nodes  plus  the  number  of 
information  nodes. 

For  some  systems  the  number  of  equations  to  be  solved  may 
be  very  large.  However,  Doerry  points  out  in  the  example 
above  that  the  system  may  be  broken  into  a  set  of  blocks  which 
may  be  solved  individually  resulting  in  a  faster  solution. 
The  above  equations  may  be  broken  into  two  1X1  blocks,  one  2X2 
block,  and  2  more  1X1  blocks  [l:p.24] .  WAVESIM's  method  of 
identifying  blocks  to  be  sequentially  solved  uses  the  system 
structural  jacobian. 

The  first  step  in  developing  a  system  structural  jacobian 
is  to  order  the  system  variables  and  system  equations.  The 
system  structural  jacobian  is  then  produced  by  combining 
device  structural  jacobians  according  to  the  arithmetic  of 
structural  jacobian  elements.  Recall  Table  2  which  lists  the 
codes  used  in  structural  jacobians.  The  arithmetic  used  to 
combine  these  elements  is  dependent  on  the  representation  of 
the  variables  of  interest  according  to  the  following 
arithmetic  rules  [l:p.58] : 

1+1=1) 

1-1=0 

I±0=I 

-J±0=D  <38> 

-I±I=D 

±n±m=±m±n=n  (nam,ni*Ij 
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The  resulting  system  structural  jacobian  now  shows  the  nature 
of  the  dependence  of  system  equations  to  each  of  the  system 
variables. 

The  system  structural  jacobian  facilitates  reducing  the 
system  of  equations  into  smaller,  more  easily  solved  blocks. 
Blocks  are  identified  in  order  and  each  block  depends  on 
system  variables  either  previously  solved  for  or  system 
variables  solved  within  the  current  block.  Doerry  provides 
algorithms  for  ordering  equations,  building  the  system 
structural  jacobian  and  reducing  the  system  into  blocks.  The 
best  reduction  results  in  the  largest  number  of  small  blocks. 
[l:pp. 59-60] 

Once  the  reduced  system  of  equations  is  built  the 
individual  blocks  are  solved.  Each  block  is  itself  a  small 
system  of  algebraic  equations  which  may  be  solved  using  an 
iterative  technique  such  as  Gauss - Seidel ,  Gauss- Jacobi,  or 
Newton- Raphson. 

3 .  Waveform  Representation  of  Interface  Variables 

The  simulation  of  dynamic  systems  using  a  digital 
computer  requires  numerical  integration  or  differentiation. 
Numerical  techniques  require  close  attention  to  the  time  step 
used  due  to  requirements  on  accuracy  and  the  potential  problem 
of  numerical  stability.  Problems  arise  if  the  system  of 
equations  is  stiff.  A  stiff  set  of  equations  is  one  with 
widely  varying  time  constants.  Unfortunately  conditions 
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unique  to  the  relatively  small  power  plant  associated  with 
shipboard  power  systems  result  in  a  tightly  coupled  the  system 
described  by  a  stiff  set  of  equations. 

WAVESIM  addresses  the  problem  of  tightly  coupled 
systems  using  waveform  representation  of  system  variables. 
The  simulation  is  broken  into  distinct  time  intervals. 
WAVESIM  adjusts  these  time  intervals  from  step  to  step 
depending  on  the  behavior  of  the  system  being  studied.  Over 
each  time  interval  the  system  variables  are  represented  as 
vectors  of  coefficients.  To  arrive  at  the  time  domain  values 
of  the  variables,  the  vector  components  are  interpreted  as 
coefficients  of  a  particular  type  of  series.  The  minimum 
information  necessary  to  convert  the  vector  of  coefficients  to 
a  time  domain  solution  is  [1:  p.63] : 

•  beginning  and  ending  times  of  interval 

•  number  of  coefficients 

•  type  of  series  used  in  the  representation 

•  a  unique  name  or  identifier 

•  the  vector  of  coefficients. 

WAVESIM  has  available  the  following  series  available  to 
represent  system  variables: 

•  Legendre  series 

•  Data  series 

•  Fourier  series 

•  Chebyshev  series 

•  Polynomial  series 
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WAVESIM  uses  the  waveform  content  of  the  highest  order 
term  in  the  series  as  a  measure  of  the  accuracy  of  the 
representation.  Waveform  content  is  simply  the  magnitude  of 
the  highest  order  term  divided  by  the  square  root  of  the  sum 
of  the  squares  of  all  the  coefficients. 

The  advantages  of  waveform  representation  over 
conventional  methods  are  [l:p.62]: 

•  Interpolation  between  time  increments  is  not  required. 
The  value  of  any  variable  may  be  readily  determined  at  any 
time. 

•  Numerical  stability  of  integration  and  differentiation  is 
not  an  issue.  Time  step  control  is  strictly  a  matter  of 
numerical  accuracy.  Integration  and  differentiation  are 
linear  matrix  operators  when  using  waveform  methods. 

•  Conversion  between  waveform  types  is  usually  a  linear 
matrix  operation  allowing  the  most  efficient  waveform  type 
to  be  used  depending  on  the  operation  being  performed. 

A  major  difficulty  arising  with  the  use  of  waveform 
operators  is  that  the  accuracy  of  waveform  representation 
decreases  sharply  if  the  function  represented  is 
discontinuous.  In  order  for  the  function  to  be  accurately 
represented,  waveform  boundaries  should  fall  at  the  time  of 
the  discontinuity.  In  simulating  a  system  in  which  the  time 
of  a  discontinuity  is  known  this  problem  is  simply  addressed 
by  inserting  a  breakpoint  in  the  input  file.  An  example  of 
this  situation  might  be  when  conducting  a  system  fault  study 
in  which  a  particular  phase  or  phases  are  grounded  at  a  known 
time.  Breakpoint  insertion  is  not  viable  if  a  system  modelled 
includes  solid  state  power  converters.  Doerry  identifies  the 
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need  for  robust  and  accurate  methods  to  predict 
discontinuities  [l:p.l63]. 

4.  Newton -Raphson  Solution  to  Systems  of  Algebraic 

Equations 

The  systems  of  possibly  non-linear  equations  resulting 
from  the  reduction  of  the  network  being  studied  may  be  solved 
using  an  iterative  technique  to  generate  a  consistent  set  of 
import  variables.  An  initial  estimate  is  made  for  the 
solution  and  the  equations  are  solved  and  compared  to  the 
known  output  resulting  in  an  error.  Based  on  this  error, 
corrections  are  generated  and  added  to  the  assumed  solution. 
This  process  is  repeated  with  the  new  assumed  conditions  until 
the  error  for  all  system  variables  is  less  than  a  specified 
threshold.  The  method  used  to  generate  the  correction  to  the 
assumed  solution  differentiates  between  the  various  solution 
methods.  WAVESIM  uses  the  Newton -Raphson  algorithm  modified 
to  include  a  continuation  parameter  as  discussed  below.  The 
basic  Newton- Raphson  method  begins  with  a  set  of  algebraic 
equations 


fy  (x) 
f2  (x) 


fix)  H  * 


=y 


fN(x) 


(38) 
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where  y  and  x  are  vectors  and  f  (x)  is  a  vector  of  functions. 
Equation  (38)  may  be  rewritten 

y-f  (x)  =0  (39) 


Making  an  initial  guess  for  x  (the  inputs)  the  corrected 
solution  is 

x(i+l)  =x(i) +J-1  (i)  [y-f  (x(i)  )  ]  (4°) 


The  matrix  J'1  used  to  correct  each  guess  is  the  inverse  of 
the  jacobian  matrix.  This  matrix  is  based  on  the  Taylor 
series  expansion  of  f (x)  about  an  operating  point  x(i) 


y=f(x(i))+~  |„(i)  (x-x(i))+JL 


cbc 


2! 


d2f 

dx2 


!x=x(i)  (x-x(i) )  2  + 


(41) 


Neglecting  higher  order  terms  and  solving  for  x 


x(i+l)  =x(i)  + 


df  | 


-1 


[y-f(x(i)J 


(42) 


By  comparing  (40)  and  (42)  J'1  is  given  by  (43) . 

One  can  observe  that  the  Jacobian  matrix  must  be  non- singular 
if  Newton -Raphson  is  to  be  used.  [20] 

There  are  multiple  reasons  why  Newton -Raphson  may  fail 
to  arrive  at  a  correct  solution  for  a  nonlinear  algebraic 


69 


■  3f,  Bfy 
dx^  dx2 
df2  df2 
dx\  dx2 


afii-1 

df2 


(43) 


0fN  dfH 
dx1  dx2 


dfH 


lx«x(i ) 


equation.  The  path  traced  by  x  for  successive  iterations  is 
referred  to  as  the  trajectory.  Around  each  solution  (there 
may  be  multiple  solutions  to  a  nonlinear  equations)  there  is 
a  region  of  convergence.  If  the  initial  guess  is  within  the 
correct  region  of  convergence  the  trajectory  will  converge  to 
a  correct  solution.  If  the  initial  guess  is  not  within  this 
region  of  convergence  one  of  five  outcomes  is  possible 
[1 :p. 72]  : 


•  The  trajectory  may  diverge 

•  The  solution,  may  by  chance,  enter  the  region  of 
convergence  hence  arriving  at  a  correct  solution 

•  The  Jacobian  matrix  may  become  singular 

•  The  trajectory  may  become  cyclic 

•  The  trajectory  could  enter  a  chaotic  region  from  which  it 
does  not  emerge,  become  cyclic,  or  converge  to  the  correct 
solution. 


The  nature  of  the  nonlinear  equations  determines  the  size  of 
the  region  of  convergence.  For  linear  equations,  the  region 
of  convergence  is  infinite. 
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WAVESIM  addresses  convergence  issues  using  a  method 
known  as  homotopy  in  which  the  continuation  parameter 
mentioned  in  the  discussion  of  device  modelling  is  used.  This 
method  attempts  to  drive  the  solution  trajectory  into  the 
region  of  convergence  in  the  following  manner.  Given  a 
nonlinear  equation 

F(x,y)  =  0  (44) 


a  new  equation  is  generated 

F'{x,y,  a)  =0 


(45) 


where 

F'(x,y,  0)  =G(x,y) 
F'{x,y,  1)  =F(x,y) 


(46) 


G(x#y)  is  a  linear  function  of  x.  WAVESIM  accomplishes  this 
by  setting 

H{x,y,  a)  =a H{x,y)  -  (1-a)  G{x,y)  (47) 

Each  nonlinear  device  model  should  include  a  linearized  set  of 
constitutive  equations  to  support  this  function.  [l:p.75] 

Through  the  use  ut  the  default  command  the  user  may 
set  the  initial  value,  the  initial  increment,  and  the  minimum 
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increment  of  a.  Two  approaches  for  using  a  to  arrive  at  a 
correct  solution  are  presented  [l:p.75] . 

One  method  begins  with  a=i.  if  the  algorithm  fails  to 
converge  to  a  solution  after  a  predetermined  maximum  number  of 
iterations,  a  is  decremented.  This  process  is  continued  until 
the  region  of  convergence  (which  is  assumed  to  increase  in 
size  with  decreasing  a)  encompasses  the  initial  guecs.  The 
parameter  a  is  then  incremented  and  theoretically  the 
trajectory  converges  on  the  correct  solution. 

The  other  method  given  in  [1]  is  to  start  with  a=0. 
Since  U(x,y)  is  a  linear  function  the  region  of  convergence  is 
infinite.  The  parameter  a  is  then  incremented  and  the  initial 
guess  is  stepped  towards  the  solution.  Doerry  demonstrates 
however,  that  this  method  is  not  always  successful  due  to 
"bifurcations  of  solutions  as  a  is  incremented."  [l:p.75] 

C.  WAVE SIM  ALGORITHM 

Figure  21  shows  the  simulation  flow  chart  for  WAVESIM 
[1 :p. 77] .  The  first  portion  of  the  program  reads  the  input 
file,  builds  the  structural  jacobian,  initializes  the 
simulation  parameters,  detects  the  block  sequence  and  reduces 
the  system.  The  next  portion  of  the  program  sequentially 
solves  the  given  blocks.  Figure  22  expands  the  "solve  blocks" 
portion  of  Fig.  21  to  show  the  steps  required  to  solve 
individual  blocks  [l:p.84].  Note  that  "success  returned"  is 
reached  only  if: 
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Figure  21  Flow  chart  of  WAVESIM  simulation  algorithm 
[1 :p. 77] 
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Figure  22  Flow  chart  for  "solve  blocks"  portion  of  Fig.  2 
[  1 : p . 84 ] 


•  the  calculated  errors  (for  the  assumed  solution)  are  small 

•  a  =  1 

•  The  truncation  error  is  within  limits 
Failure  is  returned  if: 

•  after  an  apparently  successful  convergence  the  truncation 
error  is  too  large 

•  the  solution  is  diverging  between  iterations 

•  the  allowable  number  of  iterations  is  exceeded  and  a  is 
too  small 

If  the  blocks  successfully  converge  then  the  solution  for 
the  given  time  interval  is  saved  and  the  new  time  interval  is 
computed  and  solved.  This  process  continues  until  the 
simulation  is  complete. 

If  a  block  fails  to  converge  either  the  ending  time  of  the 
current  interval  is  reduced  (hence  the  interval  duration  is 
shortened)  or  the  number  of  coefficients  is  increased  (the 
program  decides  which  action  to  take) .  If  the  specified  limits 
of  the  time  interval  control  and  number  of  coefficients  are 
not  exceeded,  the  program  will  attempt  to  solve  the  blocks  for 
the  new  conditions.  However,  if  the  simulation  limits  are 
exceeded,  the  program  terminates  and  prints  an  error  message. 
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VI  ANALYSIS  AND  TESTING  OF  WAVE SIM 


The  purpose  of  this  thesis  is  to  develop  a  set  of  metrics 
to  evaluate  software  used  to  conduct  detailed  modelling  of 
shipboard  power  systems  using  digital  computers.  Chapter  II 
section  D.3  introduces  two  efforts  in  this  area.  Chapter  III 
enumerates  a  set  of  metrics  as  derived  from  work  in  the  area 
of  detailed  digital  modelling  of  power  systems.  This  chapter 
will  evaluate  WAVESIM  in  terms  of  the  metrics  presented  in 
Chapter  III. 

A.  STUDY  APPROACH 

The  analysis  of  WAVESIM  will  be  conducted  in  two  steps. 
The  first  step  is  to  run  WAVESIM.  The  second  step  is  to 
evaluate  WAVESIM  in  terms  of  the  metrics  enumerated  in  Chapter 
III.  To  accomplish  this  the  turbine  emulator  of  the  MSL, 
introduced  in  Chapter  II  section  D.l  and  described  in  detail 
in  Chapter  IV  section  A,  will  be  modelled  using  devices 
supplied  by  LT  John  Amy,  USN,  of  the  Massachusetts  Institute 
of  Technology.  Each  component  of  the  turbine  emulator  will  be 
tested  individually  and  then  the  entire  turbine  emulator  model 
will  be  tested.  Also  the  voltage  regulator  model  with  its 
non-linear  saturation  effects  will  be  tested  to  study  the 
ability  of  WAVESIM  to  model  non-linear  devices. 
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The  response  of  the  turbine  emulator  and  its  individual 
components  to  a  pulse  disturbance  will  be  simulated. S t e  a  dy 
state  operation  with  zero  error  will  be  established  followed 
by  the  application  of  a  pulse  disturbance.  The  duration  of  the 
disturbance  will  be  adequate  to  once  again  establish  steady 
state.  Finally,  the  system  will  be  allowed  to  settle  to  steady 
state  following  termination  of  the  disturbance.  Similar 
simulations  will  be  conducted  using  the  software  package 
SIMULAB  by  The  Math  Works,  Inc.  and  the  results  will  be 
compared. 

For  both  parts  of  the  evaluation  of  WAVESIM,  modelling  the 
turbine  emulator  and  comparing  WAVESIM  to  the  metrics 
presented,  the  program  and  models  are  treated  as  a  commercial 
package.  No  effort  is  made  to  correct  problems  within  device 
models  nor  are  additional  device  objects  developed.  While  it 
is  noted  that  system  developers  would,  in  most  cases,  build 
some  of  their  own  models  to  account  for  the  uniqueness  of  a 
system  under  study,  at  this  time  WAVESIM  is  not  adequately 
developed  to  make  this  a  worthwhile  effort. 

B.  MODELING  THE  TURBINE  EMULATOR 
1 .  Speed  Regulator 

The  block  diagram  for  the  speed  regulator  is  given  in 
Fig.  9  of  Chapter  IV.  The  parameters  for  this  device  as  given 
by  Mayer  are  [14]  : 
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•  Ks=5 . 3  x  10 

•  Ks1=4 

•  rs1=.561  seconds 

•  ts2=28  x  10'3  seconds 

•  ts3=6.028  x  10"3  seconds 

•  rs 4=.9  seconds 

•  tS5=  .  45  seconds 

•  ts6=5  x  10’3  seconds 

•  Ci  =4.107  X  10‘3  Volts/Ampere 


The  test  was  conducted  with  ttref=188.5  radians/second.  The 
motor  speed  was  input  according  to: 


Cd 


188 .5 rad/sec 
94.25rad/sec 
I88.5racf/sec 


t<4  seconds 
4<  t<8  seconds 
8  <  t  seconds 


(48) 


A  fixed  time  step  of  1  ms  was  used  for  the  WAVESIM  simulation 
and  for  the  SIMULAB  simulation.  Figure  23  is  a  plot  of  the 
WAVESIM  results  superimposed  on  the  SIMULAB  results.  The 
reader  will  observe  the  discrepancies  between  the  SIMULAB  and 
the  WAVESIM  results.  WAVESIM  fails  to  accurately  depict  the 
dynamic  response  and  omits  the  fast  transients  associated  with 
the  step  inputs  resulting  in  potentially  erroneous  values  of 
the  maximum  and  minimum  currents  in  the  circuit.  Table  III 
summarizes  the  data  provided  by  Fig.  23. 
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Figure  23  Speed  regulator  response  to  a  pulse  disturbance 
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time  in  seconds 


|  TABLE  III  SPEED  REGULATOR  RESPONSE  TO  A  PULSE  DISTURBANCE 

Initial 

steady 

state 

Final 

steady 

state 

Maximum 

Iref 

Minimum 

■^ref 

WAVESIM 

-411  A 

1528  A 

1740  A 

-411  A 

SIMULAB 

-433  A 

1513  A 

2536  A 

-1197  A 

At  this  point,  an  analytic  verification  of  the  results  is 


required  to  determine  whether  SIMULAB  or  WAVESIM  is  in  error. 
Two  cases  are  derived: 

•  wm=“ref 

For  the  block  diagram  of  the  speed  regulator  given  in  Fig.  9 
the  analytic  solution  may  be  determined  by  taking  the  inverse 
Laplace  transform  of  the  given  transfer  functions  using 
partial  fraction  expansion.  For  the  case  of  Mm=o>ref=188 . 5  the 
analytic  solution  in  the  time  domain  is: 

izgf{  t)  =24344e'200t-2l272e~165t-104e'2’22t-433  (49) 

Figure  24  shows  the  simulation  results  of  WAVESIM  and  SIMULAB 
superimposed  with  the  analytic  result  for  the  conditions  given 
for  (48)  .  Note  that  SIMULAB  agrees  exactly  with  the  analytic 
solution.  The  analytic  expression  for  iref(t)  when  «ref=188.5 
rad/sec  and  um=94.25  rad/sec  is  given  by: 
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Figure  24  Comparison  of  solution  methods  for  the  speed 
regulator  response  to  a  step  input  with  «ref=«^ 
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seconds 


iI9f{  t)  =l2100e-2OOt-9600e'1S5t-52e‘a*a2t+486 .5t+67  (50) 

Figure  25  shows  plots  of  (49) ,  WAVESIM,  and  SIMULAB  results 
under  the  same  conditions.  Once  again  the  SIMULAB  output 
essentially  matches  the  analytic  solution  while  WAVESIM  yields 
inaccurate  results  for  the  transient  response. 

The  analytic  results  expressed  in  (48)  and  (49)  do  not 
include  the  amplifier  limiter  shown  in  Fig.  9.  The  SIMULAB 
models  used  include  this  feature.  One  may  infer  that  since 
the  SIMULAB  results  exactly  match  the  analytic  results,  the 
limiter  does  not  affect  these  particular  simulations. 

2 .  Current  Regulator 

The  block  diagram  for  the  current  regulator  is  given 
in  Fig.  10  of  Chapter  IV.  The  parameters  for  this  device  as 
given  by  Mayer  are  [14] : 

•  C1=4.107xl0'3  V/A 

•  Kj  =  0 . 452 

•  Kji=2 . 159 

•  Kj2=71 . 45 

•  Kj3=5 . 076xl0"5 

•  Tn=6.8xlOP'3  seconds 

•  rj2=l .  8xl0'3  seconds 

•  Tj3=7.3xl0'3  seconds 

The  simulation  is  conducted  with  iref=1000  Amperes.  The  motor 
current  is  given  by: 
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Figure  25  Comparison  of  solutions  for  speed  regulator 
response  to  a  step  input  with  wref,*«*m 
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seconds 


1000  amperes 
800  amperes 
1000  amperes 


t<  .2  seconds 
.  2<  t< . 4  seconds 
.  4  <  t  seconds 


(51) 


The  minimum  time  increment  allowed  is  1  /xsec  and  the  maximum 
is  500  nsec.  Figure  26  is  a  plot  of  the  WAVESIM  and  SIMULAB 
results.  The  WAVESIM  solution  completely  misses  the 
transients  caused  by  the  step  changes  in  input.  Further,  the 
steady  state  values  of  the  terminal  voltages  disagree  by  7.9 
volts.  Table  IV  highlights  the  information  shown  in  Fig.  26. 


TABLE  IV  RESPONSE  OF  THE  CURRENT  REGULATOR  TO  A  PULSE 
DISTURBANCE 

Initial 

steady 

state 

Final 

steady 

state 

Maximum 

Vm 

Minimum 

Vm 

WAVESIM 

0  v 

11.45  v 

11.45  v 

0  v 

SIMULAB 

-7.83  v 

3.62  v 

5.47  v 

-7.83  v 

Analytic  solutions  for  the  current  regulator  response  to 
two  input  conditions  are  given  in  (50)  and  (51)  .  For  the  case 
of  im=iref =1000  amperes  the  current  regulator  response  is 


vm{  t)  =7  .  ie~555-6t+  ,71e~137,0t-7 .83  (52) 


For  the  case  of  iref=1000  amperes  and  im=800  amperes  the 
current  regulator  response  is 

vm(t)  =5.4e-5S5-6t+.57e'137-0t+57 .3t-6  (53) 
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Figure  26  Current  regulator  response  to 
disturbance 


pulse 
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Figure  27  Comparison  of  results  for  current  regulator 
response  to  a  step  input  with  iref=im 
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Figures  27  and  28  are  plots  of  (52)  and  (53)  respectively 
superimposed  on  the  SIMULAB  and  WAVES IM  results  for  the  same 
conditions.  Once  again  there  is  exact  agreement  between 
SIMULAB  and  the  analytic  result  and  a  lack  of  agreement 
between  WAVESIM  and  the  analytic  result. 

3 .  DC  Motor 

The  DC  motor  of  the  turbine  emulator  is  rated  at  800 
hp,  500  volts  and  1750  rpm.  The  block  diagram  for  the  motor  is 
shown  in  Fig.  11  of  Chapter  IV.  The  parameters  for  this 
machine  as  given  by  Mayer  are  [21] : 

•  ra=4 . 18xl0'3  Q 

•  La= .  Ill  Mh 

•  Ta=La/ra=.266  ms 

•  Kv=2 .487  V*s/rad 

•  Ki=l . 769  f t*lb/A  =2.4  N*m/A 

The  analytic  expression  for  this  model  is 

i„<  (54) 

■*" a 


where 


err=vref-Kv*oa 


(55) 
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Figure  29  is  a  plot  of  the  SIMULAB,  WAVESIM  and  analytic 
solutions  for  w„,=188.5  and  vm=450  volts.  All  three  methods 


agree . 

4.  Combined  Elements  of  the  Turbine  Emulator 

Figure  30  shows  the  simulation  results  for  the  overall 
turbine  emulator  simulation.  For  this  simulation  the  speed 
regulator,  the  current  regulator  and  the  DC  motor  are 
connected  according  to  Fig.  8.  The  dynamic  response  of  the 
motor  current  was  simulated  with  wref=188.5  and  as  follows: 


188.5 

94.25 

188.5 


t< .  8  seconds 
. 8 < t<i . 6  seconds 
1 . 6  <  t  seconds 


The  nature  of  the  response  given  by  WAVESIM  is  completely 
different  than  the  response  given  by  SIMULAB.  Table  V 
highlights  the  data  presented  in  Fig.  30. 


TABLE  V  RESPONSE  OF  THE  TURBINE  EMULATOR  TO  A  PULSE 
DISTURBANCE 

Initial 

steady 

state 

Final 

steady 

state 

Maximum 

motor 

current 

Minimum 

motor 

current 

WAVESIM 

-343.9  A 

.72  A 

37798  A 

-87837  A 

SIMULAB 

-453.7  A 

-54.1  A 

15945  A 

-35337  A 

An  interesting  problem  occurred  when  executing  this 


simulation.  The  initial  attempt  to  conduct  the  simulation  set 
the  maximum  time  increment  to  100  ms  and  the  minimum  time 
increment  to  1  ms.  With  these  time  controls  the  program 
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Figure  30  Turbine  emulator  response  to  a  pulse  disturbance 
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failed  to  converge  to  a  solution  and  decremented  the 
continuation  parameter  to  zero.  At  this  point  the  simulation 
entered  an  apparently  infinite  loop  (ran  without  progress  for 
24  hours)  returning  the  converge  failure  warning.  The  problem 
was  avoided  by  changing  the  maximum  allowable  time  step  to  10 
ms  which  resulted  in  successful  convergence  to  a  solution 
(albeit  a  questionable  solution) . 

The  flow  charts  given  in  Figs.  21  and  22  show  that,  if  a 
block  fails  to  converge  and  a  is  too  small,  either  the  ending 
value  of  the  time  interval  should  be  updated  such  that  the 
time  interval  length  is  reduced  or  the  number  of  coefficients 
is  increased  in  order  to  decrease  the  truncation  error.  This 
did  not  happen  in  the  above  instance.  Once  alpha  was 
decremented  to  zero,  the  program  continued  attempting  to 
converge  changing  neither  the  time  interval  nor  the  number  of 
coefficients  apparently  never  leaving  the  "solve  blocks" 
portion  of  Fig.  21. 

5.  Voltage  Regulator 

The  block  diagram  for  a  voltage  regulator  is  given  in 
Fig.  12  of  Chapter  VI.  The  parameters  for  this  device  as 
given  by  Mayer  are  [14] : 

•  ke=i 

e  ka=200 

e  Kp=0 . 3 

•  ta=0.02  seconds 
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•  tf1=0 . 15  seconds 

•  rF2= 0.04  seconds 

•  A=0. 09826 

•  B=0. 33876 

The  test  was  conducted  with  vref=l  per  unit.  The  terminal 
voltage  was  input  as 

.8  pu  tcio  seconds 

V  — 

t  [1.2  pu  t>i0  seconds 

Figure  31  shows  the  results  of  this  simulation.  The  voltage 
regulator  begins  to  saturate  at  about  3.1  volts  and  is  fully 
saturated  at  3.3  volts.  There  is  general  agreement  between 
SIMULAB  and  WAVESIM  during  the  initial  phase  of  the  simulation 
with  the  two  agreeing  exactly  on  the  voltage  at  saturation. 
However,  when  the  terminal  voltage  is  stepped  to  1.2  per  unit 
there  is  delay  in  the  response  predicted  by  WAVESIM  of  over 
three  seconds. 

6.  Conclusions  on  Modelling  the  Turbine  Emulator 

The  preceding  sections  documented  efforts  to  use 
WAVESIM  to  model  the  turbine  emulator  as  presented  in  Chapter 
IV.  The  results  were  not  encouraging.  Every  effort  was  made, 
short  of  disassembling  the  models  themselves,  to  achieve 
results  that  could  be  corroborated  by  SIMULAB.  The  SIMULAB 
models  were  verified  against  analytic  solutions  to  ensure 
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Figure  31  Voltage  regulator  response  to  a  step  change  i 
terminal  voltage 
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time  in  seconds 


their  accuracy.  In  every  case,  WAVESIM  failed  to  accurately 
reflect  the  dynamic  performance  of  the  systems  being 
simulated. 

C.  EVALUATION  OP  WAVESIM  VERSUS  THE  METRICS  IN  CHAPTER  THREE 

1.  System  and  Component  Level  Modelling  Environment 

WAVESIM  is  highly  modular.  With  properly  developed 
models  and  with  adequate  accompanying  documentation  the  user 
could  be  expected  to  quickly  develop  the  skills  necessary  to 
use  the  program.  This  assumes  that  adequate  models  have  been 
developed.  In  its  current  stage  of  development  the  available 
models  are  inadequate  for  virtually  any  task. 

Model  development  is  not  intuitive.  [21]  and  [22] 
detail  the  steps  involved  in  encoding  a  three  phase 
synchronous  generator.  The  requirement  to  cast  all 
relationships  into  waveform  equations  is  burdensome.  The  use 
of  waveform  equations  and  special  waveform  functions  to 
perform  even  simple  arithmetic  tasks  should  be  made 
transparent  to  the  user.  The  modelling  technique  requires  an 
intimate  knowledge  of  the  internal  workings  of  WAVESIM  making 
it  difficult  for  the  person  not  involved  in  developing  the 
program  to  write  their  own  models. 

2.  Robustness  When  Simulating  Non-Linear  or  Rapid 

Switching  Topologies 

There  are  inadequate  models  available  to  determine  if 
WAVESIM  can  model  widely  varying  topologies.  A  very  large 
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proportion  of  the  load  on  an  IED  system  will  be  solid-state 
power  converters.  Therefore,  the  ability  to  conduct  detailed 
modelling  of  solid  state  power  converters  is  essential.  Yet, 
the  methods  used  by  WAVES IM  may  be  ill  suited  for  this  kind  of 
simulation.  Doerry  notes  that  there  are  difficulties  modelling 
discontinuous  functions  and  improved  methods  of  discontinuity 
prediction  are  required  [l:p.l63]. 

3 .  Correctness  of  Solution 

WAVESIM  uses  generally  accepted  solution  methods. 
However,  as  a  modular  program  WAVESIM' s  success  is  dependent 
on  correct  development  of  component  models.  WAVESIM  is  capable 
of  arriving  at  correct  solutions  to  both  linear  and  nonlinear 
modelling  problems  as  demonstrated  in  [1]  .  Chapter  VI  of  this 
thesis  shows  that  WAVESIM  can  also  deliver  incorrect 
solutions.  A  good  deal  of  independent  validation  is  required 
on  all  models  developed  if  WAVESIM  is  to  be  trusted  as  a 
useful  simulation  tool. 

4.  Software  Domain 

Currently,  field  test  data  can  be  used  in  conjunction 
with  WAVESIM  due  to  the  programs  dependence  on  MATLAB.  While 
no  method  exists  to  use  field  test  data  to  drive  a  particular 
model,  it  should  be  possible  to  develop  an  object  to  support 
this  feature.  With  its  current  modular  construction  and  with 
its  use  of  parameters  the  ability  exists  to  vary  control 
parameters  in  support  of  control  system  synthesis.  However, 
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the  time  required  to  reach  a  new  solution  may  make  this 
process  too  time  consuming. 

5.  Implementation  Options 

WAVESIM  is  written  in  standard  ANSI  C  and  should, 
therefore,  be  highly  portable.  With  its  current  dependence  on 
MATLAB  one  may  foresee  a  potential  for  compatibility  problems 
with  future  versions  of  MATLAB. 

Long  term  maintainability  is  an  important  issue. 
WAVESIM  is  essentially  a  one  man  project.  While  the  code  is 
documented  (as  clearly  as  C  code  can  be) ,  finding  people  or 
preferably  an  organization  qualified  and  willing  to  further 
develop  and  maintain  the  program  over  the  long  term  may  be 
difficult . 

6.  Ease  of  Use 

WAVESIM  is  neither  intuitive  nor  particularly  user 
friendly.  The  file  mode  of  input  is  difficult  to  master  and 
the  large  number  of  "default"  sub  commands  as  described  in 
Chapter  Y  section  5B  makes  input  file  development  cumbersome. 
An  interactive  (preferably  graphical)  user  interface  is 
needed. 

It  is  necessary  to  specify  which  system  variables  are 
of  interest  in  the  input  file.  Only  those  system  variables 
specified  in  the  input  file  plot  command  are  available  as  data 
series  at  the  conclusion  of  the  simulation.  Specifying  more 
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system  variables  for  output  than  required  is  costly  in  terms 
of  both  simulation  time  and  memory  usage. 

7.  Software  Speed  Versus  System  Complexity 

While  only  relatively  simple  systems  were  used  for 
this  thesis,  it  is  apparent  that  execution  of  the  simulation 
file  within  MATLAB  is  slow.  The  simulations  conducted  for 
this  chapter  were  conducted  on  a  Sun  Sparc  Station  II  running 
Open  Windows  version  2  and  PROMATLAB  version  3.51.  It  was 
common  for  WAVESIM  simulations  to-  require  a  number  of  hours. 
The  SIMULAB  simulations,  with  the  same  error  tolerances  and 
identically  defined  time  steps,  took  less  than  three  minutes 
(usually  much  less) .  The  simulation  of  the  turbine  emulator 
presented  in  section  B4  of  this  chapter  required  41  minutes 
for  WAVESIM  to  complete  and  only  3  seconds  for  SIMULAB  to 
complete.  A  more  telling  comparison  is  that  WAVESIM  required 
969.6  million  floating  point  operations  (MFLOPS)  to  complete 
the  simulation  while  SIMULAB  required  only  25800  floating 
point  operations  for  the  turbine  emulator  simulation. 

8 .  Continued  Support 

As  discussed  in  section  B5  of  this  chapter,  long  term 
maintainability  is  a  primary  consideration  in  choosing  a 
software  tool.  WAVESIM  is  basically  the  result  of  a  very 
small  group  of  people  who  are  no  longer  with  the  institution 
at  which  the  program  was  developed.  While  WAVESIM  is  of 
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interesting  academic  value,  it  no  longer  has  the  support 
necessary  to  make  it  into  a  useful  tool. 
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VIII  CONCLUSIONS 


A.  NEED 

The  Navy  has  a  definite  requirement  for  a  simulation  for 
use  in  developing  and  studying  advanced  marine  power  systems. 
This  need  is  brought  about  primarily  due  to  the  pending  shift 
from  the  traditional  engineering  plant  configuration  of 
separate  propulsion  and  electrical  systems  to  an  advanced 
marine  power  system  integrating  propulsion  and  power 
generating  prime  movers. 

The  simulation  program  must  meet  the  requirements 
enumerated  in  Chapter  IV.  Specifically  it  must  be: 

•  modular  and  have  an  adequate  model  library 

•  robust  and  capable  of  simulating  rapidly  switching 
topologies  associated  with  large,  solid-state  power 
converters 

•  capable  of  providing  correct  dynamic  and  steady  state 
performance  data 

•  portable  and  maintainable 

•  easy  to  use 

•  fast  and  capable  of  solving  stiff  systems  of  equations 

•  well  supported 
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B.  ANALYSIS  OF  WAVES IM 


WAVESIM  introduces  and  combines  advanced  simulation 
techniques  in  an  effort  to  develop  a  simulation  tool 
specifically  tailored  to  modelling  advanced  marine  power 
systems.  The  result  has  been  a  test  program  that  demonstrates 
these  methods  but  in  its  current  state  is  not  useful  as  an 
actual  tool.  This  conclusion  is  based  on  the  metrics  presented 
in  Chapter  IV  and  on  the  evaluation  or  WAVESIM  as  presented  in 
Chapter  VII. 

C.  FUTURE  WORK 

1 .  WAVESIM 

Additional  testing  of  WAVESIM  including  detailed 
analysis  of  previously  developed  models,  investigating  why  the 
dynamic  performance  of  the  given  models  does  not  match  the 
analytic  results,  and  developing  new  models  are  areas  of 
future  work.  Further,  investigation  of  methods  of 
discontinuity  prediction,  improved  numerical  efficiency,  and 
removal  of  the  dependence  on  MATLAB,  as  well  as  developing  an 
online  user  interface  are  areas  of  future  work  if  Navy 
interest  remains  in  improving  WAVESIM. 

2 .  Other  programs 

The  usefulness  of  two  commercially  available  programs 
mentioned  previously  in  this  thesis  warrant  further 
investigation.  ACSL  introduced  in  Chapter  II  section  D.3.a 
has  been  used  for  both  component  and  system  level  modelling  of 
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power  systems  at  Purdue  University  [6] .  SIMULAB  was  shown  in 
Chapter  VII  to  accurately  model  the  components  of  the  turbine 
emulator.  ACSL  research  is  currently  underway  at  Purdue. 

SIMULAB  meets  many  of  the  requirements  given  in 
Chapter  III  for  a  simulation  tool.  It  is  modular  and  allows 
for  the  modelling  of  both  linear  and  non-linear  devices.  It  is 
commercially  supported.  It  is  fast  and  efficient  allowing 
various  integration  methods  to  be  specified  depending  on  the 
nature  of  the  system  being  modelled.  Lastly,  it  is  intuitive 
in  that  it  uses  a  graphical  user  interface  displaying  actual 
transfer  functions  for  the  various  blocks.  Investigation  of 
SIMULAB  to  meet,  the  need  of  a  power  systems  simulation  tool 
may  be  a  worthwhile  project. 
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